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Abstract. 

We use the Negele-Vautherin density matrix expansion to derive a quasi-local 
density functional for the description of systems of fcrmions interacting with short- 
ranged interactions composed of arbitrary finite-range central, spin-orbit, and tensor 
components. Terms that are absent in the original Negele-Vautherin approach owing 
to the angle averaging of the density matrix are fixed by employing a gauge invariance 
condition. We obtain the Kohn-Sham interaction energies in all spin-isospin channels, 
including the exchange terms, expressed as functions of the local densities and their 
derivatives up to second (next to leading) order. We illustrate the method by 
determining the coupling constants of the Skyrme functional or Skyrme force that 
correspond to the finite-range Gogny central force. The resulting self-consistent 
solutions reproduce the Gogny-force binding energies and radii within the precision 
of 1-2%. 



PACS numbers: 21.60.Jz, 21.10.Dr, 21.30.-x 



1. Introduction 

The search for a universal energy density functional (EDF) [TJ, that would be able to 
provide a spectroscopic-quality [2] description of atomic nuclei, is at the focus of the 
present-day studies in nuclear structure. Recently, several advanced phenomenological 
analyses were aimed at improving the standard relativistic [3] or nonrelativistic local 
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[U El El El E] or nonlocal functional [S] . A significant ongoing effort is also devoted to 
microscopic derivations of the functionals (see, for example, Refs. j9j El E].) One of 
the central points of the current EDF studies is the question: to what extent can the 
finite-range effective interactions be approximated as quasi-local density functionals? 

The framework to build the quasi-local theory was set up in the seminal paper 
by Negele and Vautherin (NV) [12], which introduced the so-called density matrix 
expansion (DME) method. Later, other methods to achieve the same goal, like the 
semiclassical expansion (T3J E] were also proposed and studied. The original NV 
expansion, for the scalar density, allowed for treating only the even-order terms in 
relative coordinates, and thus was applicable only to even-power gradient densities. 
In the present study, we revisit the NV expansion by adding the odd-power gradient 
densities through the local gauge- invariance condition. There are two other important 
differences of the present approach with respect to the NV original, namely, (i) we treat 
the spin (vector) densities analogously to the scalar densities and (ii) we apply a different 
DME for the direct terms. These differences are motivated by the effective-theory [15] 
interpretation of the DME advocated in the present study. The ultimate test of these 
ideas can only be obtained by analysing microscopic properties of nuclear densities [TJ] . 

We present a complete set of expressions in all spin-isospin channels, applicable to 
arbitrary finite-range central, spin-orbit, and tensor interactions. Our study is focused 
on applying the NV expansion to the Gogny interaction [161 E]- This allows us to look 
at the correct scale of the interaction range that properly characterizes nuclear low- 
energy phenomena. By performing expansion up to second order (or next-to-leading 
order, NLO), one obtains the local Skyrme functional [T8| IT9]. In this way, we establish 
a firm link between two different, local and nonlocal, very successful functionals. 

The paper is organized as follows. In Section |2l we introduce reformulation of the 
NV expansion in terms of particles without spin. Discussion of this example allows 
us to present details of the approach without complications that otherwise could have 
obscured the main ideas. In Sections El and HI we present results with spin and isospin 
degrees of freedom reintroduced, with Section H] containing applications of the formalism 



to the Gogny force. Conclusions are given in Section El and Appendix A contains 
the discussion of spin and isospin polarized nuclear matter. Preliminary results of the 
present study were published in Ref. [20J. 

2. Local energy density for spinless particles of one kind 

In this section, we consider the simplest (and academic) case of fermions with no spin 
and no isospin. First we recall that for an arbitrary non-local finite-range interaction 
V(r[, r' 2 ; r±, r 2 ), the Kohn-Sham interaction energy [21] has the form 

£int = i J d 3 r , 1 d s r' 2 d 3 r 1 d s r 2 V{r' 1 ,r' 2 ;r 1 ,r2) x 

00i, 7"'i)p( r 2, r' 2 ) - p(r 2 , r'Jpin, r 2 )), (1) 
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whereas for a local interaction, 

V(r[, r 2 ; r ls r 2 ) = ^-n^-r^K r 2 ), (2) 
the interaction energy reduces to: 

£int = i J d 3 r 1 d 3 r 2 y(r 1 ,r 2 )(p(r 1 )p(r 2 ) - p(r 2 , n)^, r 2 )), (3) 

where p(t*i) = p(ri,ri) and p(r* 2 ) = /7(t*2,t*2) are local densities. As is well known, the 
first term in Eq. (J3J) (the direct term) depends only on local densities, whereas the second 
one (the exchange term) depends on the modulus squared of the non-local density. This 
markedly different structure of the two terms requires separate treatment, as discussed 
in the following two subsections. 



2.1. Direct interaction energy 

In nuclei, the range of interaction is significantly smaller than the typical scale of the 
distance at which the local density varies. Therefore, we may expand the local densities 
p(ri) and p(r 2 ) around their average position, and use this expansion to calculate the 
direct term in Eq. 

Denoting the standard total (R) and relative (r) coordinates and derivatives as 

R=l(r 1 + r 2 ), r = r l -r 2) (4) 

V = — — — d- — -- (— —) (5) 

OR dri <9r 2 ' dr 2 \dri dr 2 ) ' 

we have the expansion of local densities to second order, 

p( ri ) = P {R + \r) = p{R) + \r a V aP {R) + |r a r 6 V Q V 6 p(H) + . . . , (6) 
p(r 2 ) = p{R - \r) = p{R) - \r a V a p(R) + \r a r h V a V b p(R) + ..., (7) 

and hence 

p(n)p(r 2 ) = p 2 (R)+lr a r b (p(R)V a V bP (R)-[V a p(R)][V b p(R)])+. . . , (8) 

where we implicitly assumed the summation over the repeated Cartesian indices a and 
b. 

Assuming that the local potential V r (r 1 ,r 2 ) depends only on the distance between 
the interacting particles, V(ri,r 2 )=V(\ri—r 2 \)=V(r), the direct interaction energy is 
given by the integral of a local energy density W^.(R), 

£g = / d s Rn^(R), (9) 
where up to second order, 



-l/int/ rf\ 1 

n dir\ JX ) — 2 



V p 2 + ±V 2 (pAp-(Vp) 2 )}+..., (10) 

where the coupling constants, V and V 2 , are given by the lowest two moments of the 
interaction, 

V n = J d 3 rr n V{r) = 4vr J drr n+2 V{r). (11) 
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After integrating by parts, the two second-order terms in Eq. fllOp are identical, so we 
can equally well use: 

Uf T {R) = \ [VoP 2 + lV 2P Ap] + . . . , (12) 

We see that the separation of scales between the range of interaction and the rate of 
change of the local density leads to a dramatic collapse of information that is transferred 
from the interaction potential to the interaction energy. Namely, the two constants, Vo 
and V 2 , completely characterize the interaction in the direct term, and the detailed form 
of the potential V(r) becomes irrelevant. Moreover, it can be easily checked that in this 
approximation, the direct energy density is exactly equal to that corresponding to the 
contact force corrected by the second-order gradient pseudopotential, 

V(r) = V S(r)-lV 2 8-6{r)d. (13) 
2.2. Exchange interaction energy 

In the exchange term of Eq. Q, the non-zero range of the interaction probes the non- 
local space dependence of the density matrix. For short-range interactions, one can 
expand p(R,r) to second order with respect to the variable r, which gives 

p(r u r 2 ) = p(R, r) = p(R) + r a d a p(R, r) + \r a r b d a d b p{R, r) + . . . 

= p(R) + ir a j a {R) + \r a r b [\V a S7 b p(R) - r ab (R)] + . . .(,14) 

where derivatives di are always calculated at r a =0, and therefore, the result can be 
expressed in terms of the standard current and kinetic densities [22J: 

j a {R) = \d a p(R,r) r=0 , r ab (R) = V^V { b 2) p(r u r 2 ) ri=r2 . (15) 

This parabolic approximation does not ensure that p(ri, r 2 ) — >0 for large r = \r\ = 
\r±— r 2 \. In the spirit of the DME [12], one can improve it by introducing three functions 
of r, vr (r), ^(r), and n 2 (r) [2U] that vanish at large r, i.e., we define the quasi-local 
approximation of the density matrix by: 

p(r u r 2 ) = n (r)p(R) + 7r 1 (r)r a d a p(R, r) + ±ir 2 (r)r a r b d a d b p(R, r) + . . . 

= n (r)p(R) + mi{r)rJap{R) + \v 2 (r)r a r b [\V a V b p{R) - r ab (R)] + . . .(16) 

Such a postulate has to be compatible with the Taylor expansion of Eq. f }T4"|) . which 
requires that 

7T (0) = 7Ti(0) = 7T 2 (0) = 1 and 7l£(0) = 7ii(0) = <(0) = 0. (17) 

Of course, for 7r (r)=7ri(r)=7r 2 (r)=l, one reverts to the parabolic approximation (fT4j) . 

The product of nonlocal densities in the exchange integral of Eq. (J3J) to second order 
reads 

pir^r^pir^r!) = nl(r)p 2 (R) 

+ n {r)7r 2 (r)r a r b {p(R)d a d bP (R, r) - [d a p{R, r)} [d b p(R, r)]} + . . . 

= 4(r)p\R) 

+ no{r)ii 2 (r)r a r b {\ P (R)V a V b p{R) - p(R)r ab (R) + Ja (R) 3b (R)} + . . . , (18) 
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where we have introduced the supplementary condition: 

vr 2 (r) = 7r (r)7r 2 (r). (19) 

This condition ensures that the quasi-local approximation of Eq. (1161) is compatible 
with the local gauge invariance [23] . Indeed, the left-hand side of Eq. (fl8|) is manifestly 
invariant with respect to the local gauge transformation, 

p(r 1 ,r 2 ) = e^-^p(r 1 ,r 2 ), (20) 

and only the difference of terms in the curly brackets in Eq. ffl8l) is invariant with respect 
to the same transformation [22j IT9] . 

Functions vr (r), vri(r), and 7r 2 (r) also depend on the parameters defining the 
approximation (Tl6|) . In particular, when the infinite matter is used to define functions 
^o{ r )i 7r i( r ) ) and ^(r), like in the DME, they parametrically depend on the Fermi 
momentum kp. By associating the local density p(R) with kp , functions 7r (r), 7r 1 (r), 
and ^(r) become dependent on p(R), and hence the damping of the density matrix in 
the non-local direction r can be different in different local points R. However, in order 
to keep the notation simple, we do not explicitly indicate this possible dependence on 
density. 

Within the quasi-local approximation, one obtains the exchange interaction energy, 
£*£ = / d 3 RH^ c (R), (21) 
where up to second order, 

KfcOR) = -HO 2 + K^lPAp - (pr - j 2 ))] + . . . , (22) 

and where r = r aa . The coupling constants, V®q and V®$, are given by the following 
moments of the interaction, 

Vg = f d 3 rr n n(r)Trj{r)V(r) = 4vr / drr n+2 7r i (r) 7 r i (r)F(r). (23) 



Unlike the coupling constants defining the direct term (fTTj) . these in Eq. (1231) have to be 
understood as the running coupling constants, which in the DME depend on the scale 
of the Fermi momentum kp or density p(R). 

Again, the separation of scales between the range of interaction and the rate of 
change of the density matrix in the non-local direction results in the dependence of 
the local energy density on two coupling constants only, and not on the details of 
the interaction. For the parabolic approximation of Eq. (THj) . the coupling constants 
that define the direct and exchange energies are identical, i.e., V®q=Vq and V^ 2 2= ^2; 
however, for the quasi- local approximation of Eq. (fTEj) they are different. This important 
observation is discussed in Sec. 13.21 in more detail. 

In nuclei, the separation of scales discussed above is not very well pronounced. 
The characteristic parameter, which defines these relative scales, is equal to kpa, where 
a stands for the range of the interaction. For example, within the Gogny interaction, 
which we analyze in detail in Sec.HJ there are two components with the ranges of a = 0.7 
and 1.2 fm, whereas kp = 1.35 fm -1 , which gives values of kpa = 0.95 and 1.62 that are 
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dangerously close to 1. Therefore, in the expansion of the local energy density ( 12 2 p . one 
cannot really count on the moments V* J n decreasing with the increasing order n. 

Instead, as demonstrated by Negele and Vautherin [12], one can hope for 
tremendously improving the convergence by using at each order the proper counter- 
terms, which make each order vanish in the infinite matter. Without repeating the 
original NV construction, here we only note that the net result consists of adding and 
subtracting in Eq. (fl"6l) the infinite-matter term, that is, 

p(ri,r 2 ) = u (r)p(R) + iu 1 (r)r a j a p(R) 

+ \v2{r)r a r b [\V a V b p{R) - r ab (R) + \5 ab k 2 F p{R)\ +... , (24) 

where we have defined functions z/j(r) such that: 

M r ) = Ko{r) - jQ(k F r) 2 7r 2 , z/i(r) = 7Ti(r), v 2 {r) = vr 2 (r). (25) 

By neglecting the term quadratic in u 2 , we can now use the approximation ( 124"1) to 
calculate the product of densities in Eq. (Tl8|) . which gives 

p^r^p^, n) = v 2 (r)p 2 (R) + u (r)u 2 (r)r a r b {\p(R)V a V b p(R) - p(R)r ab (R) 

+ \5 ab k 2 F p 2 (R) + ]a (R) 3b (R)) + . . . , (26) 

where again the gauge invariance requires that 

ul{r) = u Q {r)u 2 {r). (27) 

We note that the gauge- invariance conditions for the functions 7Tj(r) and z^(r), Eqs. (f!9l) 
and (T2~Tj) . are compatible with one another only up to the term vr|(r), which was shifted 
to higher orders. Finally, approximation ( l26l) gives the energy density analogous to ( l22l) . 

= -|[K°oV + |K° 2 2 (ipAp - ( P r-j 2 ) + lk 2 F p 2 )} +..., (28) 

where the coupling constants, V^ and V® 2 2 , are given by the moments of the interaction 
calculated for functions ^i(r), namely, 

Vll = [ dVrXrji/jM^fr) = 4vr f dr r^v^v^Vir). (29) 



We see that the two sets of auxiliary functions, Hi(r) and z/j(r), are suitable for 
discussing the approximate forms of the nonlocal density p(r* 1; r 2 ) and exchange energy 
density W^ C (R), respectively. Although for the nonlocal density they correspond to 
a simple reshuffling of terms, which gives relations (|25|) between 7Tj(r) and z^(r), for 
the exchange energy density they constitute entirely different approximations, given in 
Eqs. (1221) and ( 1281) . with expansion ( |28l) having a larger potential for faster convergence. 
Only this latter expansion is further discussed. 



2.3. Determination of functions vrj(r) 

For a given nonlocal density p(ri, r 2 ), the auxiliary functions 7r.;(r) or z/i(r), which define 
its quasi-local approximation, can be calculated as their best possible approximations 
in terms of local densities. However, the usefulness of the expansion relies on the 
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assumption that generic forms of these functions can be estimated and then applied 
to all many-body systems of a given kind. 

The standard Slater approximation j2H [25] , which is routinely used to evaluate the 
Coulomb exchange energy (cf. Refs. [261 [27]), corresponds to 

vro(r) = u (r) = and vr 2 (r) = 0. (30) 

The NV expansion [T2] gives a second-order estimate by making the momentum 
expansion around the Fermi momentum k F of an infinite system. This gives: 

= MM-) + 2i j3 (^) M (31) 

ov ; 2k F r 504 ' v ; 

105j 3 (**r) {k F rf {k F rf 

Mr) - {k F rf ~ l 18 + 792 ' {62} 

0V ; k F r 10 280 ' K J 

where j n (k F r) are the spherical Bessel functions. 

The NV functions vr (r), 7r 2 (r), and z/ (r) are plotted in Fig. [I] with solid, dashed, 
and dotted lines, respectively. One can see that for large k F r, functions 7r (r) and 7r 2 (r) 
have zeros close to one another and the same signs. Indeed, asymptotically both behave 
like cos(k F r)/k F r. Therefore, the gauge-invariance condition (j!9p can be satisfied almost 
everywhere. On the other hand, functions ^o(r) and z^ 2 (r) = ^(r), have asymptotically 
opposite signs, and the corresponding gauge-invariance condition ( )27|) can be almost 
nowhere satisfied. Nevertheless, as discussed in Sec. I2.2[ what really matters are the 
moments of interaction (1231) and (|29"1) . where functions VTj(r) and z^(r) are probed only 
within the range of the interaction, that is, up to k F r ~ 1-2. Therefore, for the NV 
expansion, one can safely use approximations: 

ni (r)= + y/\Mr)*2(r)\, (34) 



Mr) = +^Wo(r)Mr)\, (35) 

which are valid up to the first zero of j'3 or j\, respectively, that is, up to k F r ~ 7.0 and 
4.5. 

By the same token, we can replace in Eqs. (l3~T i) -( l33l the Bessel functions by 
Gaussians having the same leading-order dependence on k F r, namely, 

Mr) = exp (-^) , (36) 

?r 2 (r) = exp — — — , (37) 




P(-^]----p(-^-^). (33) 

As seen in the bottom panel of Fig. (TJ in this way, in the region of small k F r, one obtains 
a very good reproduction of the NV functions 7Tj(r) and z/j(r). 
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Figure 1. Dependence of the functions 7ro(r) (solid lines), iT2{r) (dashed lines), and 
i"o(r) (dotted lines) on kpr. The top and middle panels show the NV functions 
of Eqs. (|3T |) -([33 | . with the top panel plotted in expanded scale to better show the 
details at large ftp r. The bottom panel shows the functions 7To(r), ^(r), and vo{r) 
approximated by Gaussians as in Eqs. (|36 |) -([38 |) . 



3. Local energy density for particles with spin and isospin 

3.1. Density matrix expansion with spin and isospin 

For nucleons, the density matrix p(r\0~iTi,r 2 o~ 2 T 2 ) depends not only on positions r\ 
and r 2 but also on spin o"i,o"2 = ±1 and isospin ti,t 2 = ±1 coordinates. Since the 
strong two-body interaction is assumed to be isospin and rotationally invariant, it is 
convenient to represent the standard density matrix p{riaiTi,r20~2T2) through nonlocal 
densities p^ifi, ^2) as: 

3 

p{r 1 o- 1 T 1 ,r 2 o-2T 2 ) = j P^k(ri,r 2 )(a 1 \a^\a 2 )(Ti\T k \T2), (39) 

fi=0,x,y,z k=0 
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where gq (o'x,y,z) an d to ( r i,2,3) are the unity (Pauli) matrices in the spin and isospin 
coordinates, respectively. For /i = and p, = x,y, z, the densities are scalars and 
vectors, respectively, and for k = and k = 1,2, 3, they are isoscalars and isovectors, so 
altogether the density matrix is split into the four standard spin-isospin channels. 

For the direct term, we can proceed as in Sec. 12.11 by making the Taylor expansions 
of local densities (at ri = r 2 ) in each spin-isospin channel; that is, similarly as in Eqs. ([6]) 
and ([?]), we have 

Pflk (R ± \r) = Pflk (R) ± \r a V a pA R ) + \r a r b V a V b p^ k (R) + ... . (40) 



For the exchange term, Sec. I2.2[ the analogous Taylor expansions of nonlocal densities, 
similarly as in Eq. ( fl4l) . read 



PlMk (R, ±r) = p^ k {R) ± ir a j^ ak (R) + \r a r b [\V a V b p ^(R) - T^ abk (R)\ + ... , (41) 

where the current (jfj, a k(R)) and kinetic (r^ abk (R)) densities are defined in each channel 
as in Eqs. ( IT5]) . namely, 

jfiak(R) = \d a p^ k (R,r) r=0 , T^ iabk (R) = V^ ) V ( b 2) p flk (r 1 ,r 2 ) ri =r 2 - (42) 

The local density approximation of densities in all channels, analogous to Eqs. ffTol) 
and (124"]) . is now postulated as 

p^ k (R, ±r) = 7r (r)p /lfc ( J R) ± ivri(r)r a j Mafe ( J R) 

+ \^2{r)r a r b \V a V b p^ k (R) - T^ abk (R) + ... , 



and 



Pt+kiR, ±r) = i> (r) PlJLk (R) ± ivi{r)r a j^ ak p(R) 

+ \v<i{r)r a r b |V a Vbp M fe(i?) - T^ Labk (R) + \5 ab k 2 Fpilk (R) 
At this point, we have assumed that functions 7Tj(r) and z/j(r) 



(43) 



+ . . . . (44) 
are channel- 



independent, that is, that they are scalar-isoscalar functions. In Appendix A we discuss 
this point in more detail, and we show that the postulate of simply channel-dependent 
functions 7Tj(r) and z^(r) is incompatible with properties of infinite matter, whereas the 
proper treatment of the problem leads immediately to the channel mixing and the energy 
density, which is not invariant with respect to rotational and isospin symmetries. This 
question certainly requires further study, whereas at the moment, a consistent approach 
can only be obtained by assuming the scalar-isoscalar functions 7Tj(r) and z/j(r). 

We can now apply derivations presented in Sees. 12.11 and 12.21 to the general case 
of an arbitrary finite-range local nuclear interaction composed of the standard central, 
spin-orbit, and tensor terms: 

V(n, r 2 ) = W(r) + B{r)P a - H(r)P T - M(r)P a P T 



P{r) + Q(r)P T 



L-S + 



where r=|r| 



=|ri- 

Prr 



r 2 



+ 

and 

i(l + <7l -cr 2 ), P T = §(l+T 1 of 2 
-ihrxd, S =§(01+02), S12 



R(r) + S(r)P T 



S 



12- 



(45) 

(46) 
-0V<r 2 .(47) 



The Negele- Vautherin density matrix expansion applied to the Gogny force 



10 



After straightforward but lengthy calculations, one obtains the interaction energy in the 
form of a local integral, analogous to that for the Skyrme interaction 

3 



£ 



hit 



fc=0 



C p t p\ + C t Ap p k A Pk + CU Pk r k - jl 



&k ' T k -iabk-^abk 



+ Cf(s k 



1 I I _ l I I 

2 J aak J bbk o-'abk^bak 



+ Cf J (p fe V • J k + s k ■ (V x j k )^ 

T ak — T a fef)fc, F n 



+ C? s (V-s k 



where p k = po k , r k = r bbk, jak = joak, s ak = p ak , T ak = T a bbk-, F ak — ^(rb a bk + Ubakj 
S ahk = j abk , and J ak = e a bcjcbk are the standard local densities. The isoscalar (t = 0) and 
isovector (t = 1) coupling constants C t correspond to k = and A; = 1,2,3, respectively. 

The coupling constants of the local energy density ( 1481) are related to moments of 
the interaction in the following way: 



(4f 



Co 



( 4 



Vo 



2 

2 




-2 
-2 





-1 \ 



/ W + M™ 



1 /If 02 p 



^0 + #°o + 5#°2 fc i 



^2 a F 



(49) 



96 



Co 
CI 

°o 
Cf 

°o 
Cf 

^o Vs 



•1 4 

0-1 

4 

4 
0-1 



2-4 2-2 
2-4 0-2 
8 0-80 
8 
4 2-2 



0-1 0-2 

4 0-8 

4 











4 \ 



16 


0-4 1-2 2 

1-20 

0-4 -8 

0-4 

12 24 

12 
0-12 3-6 6 

3-6 / 



f Wo \ 



B 2 

R02 

H 2 

rr02 
ti v2 

M 2 
Mf 2 

HI 



(50) 



24 



CJ J 




All the coupling constants of the local energy density 
following moments of potentials: 

X n = J d 3 rr n X{r) = 4tt J drr n+2 X{r), 



(51) 

depend linearly on the 



X*3 
un 



d rr n Vi(r)vj(r)X(r), = Air / drr n+ 2 Vi{r)vj{r)X{r) 



(52) 
(53) 
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where X stands for W, B, H, M, P, Q, P, or S. 

Again we see that whenever expansions of density matrices, Eqs. (j4*0]) and ( T44|) . 
are sufficiently accurate within the ranges of interactions, information about these 
interactions collapses to a few lowest moments. Short-range details of these interactions 
are, therefore, entirely irrelevant for low-energy characteristics of nuclear states. This 
is typical of all physical situations, where scales of interaction and observation are well 
separated, as specified in the effective field theories. The energy density characterizing 
the low-energy effects is local and depends on local densities and their derivatives up to 
second order, whereas the dynamic information is contained in a few coupling constants. 

Moreover, the detailed large-r dependence of auxiliary functions Vi(r) on position r 
is also irrelevant, because all that matters are moments (153|) which define the coupling 
constants ( )49l) -(l5Tl) describing the exchange energy, and these are influenced only by the 
small-r properties of functions fi{r). Finally, the most important feature is the kp or 
density dependence of i^(r), which determines the density dependence of the coupling 
constants. 



3.2. Local energy density corresponding to the Skyrme force 



In general, the number of moments entering Eqs. (l49l)-(l5TT) is higher than the number 
of coupling constants, and all the coupling constants are independent. However, it is 
extremely instructive to check what happens in the vacuum limit of kp = 0. This 
situation is obtained by setting z^(r) = 1, which gives the direct and exchange moments 
equal to one another, namely, X l J n =X n) and the coupling constants of Eqs. P9|) - fl5T|) 
collapse to: 
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(56) 



These coupling constants correspond exactly to those obtained for the Skyrme force (see 
Ref. [ID] for the notations and conventions used), namely, 



to = W + M 
h = -\{W 2 + M 2 ) 
*2 = \{W 2 -M 2 ) 

^ e = TE(S 2 — R 2 ) 

w = -Up 2 + q 2 ) 



toXo 

t 2 x 2 

to 



Bo + H , 
-\{B 2 + H 2 ) 
\{B 2 -H 2 ) 
je(S 2 + R 2 ) 



(57) 
(58) 
(59) 
(60) 
(61) 



The same relations are also obtained by using in the exchange term the pure 
Taylor expansions f )4T|) ; that is, by setting 7r(r) = 1, which gives X^ n =X n , and by 
using the classification of terms as in Eqs. (12"T|) - (|2"3|) . This second way of obtaining the 
approximate coupling constants leads to results independent of k F , which are, of course, 
identical to those obtained at kp = above. 

Relations (|54l)-(l56l) imply that the coupling constants of the energy functional ( 1481) 
are dependent of one another, and in fact, half of them determines the other half. This is 
exactly the situation encountered when the energy density is calculated for the Skyrme 
interaction. Then one obtains (cf. Ref. [23J): 
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C VJ = 3C? J . (64) 

It is obvious that the above relations among the coupling constants result from an 
oversimplified approximation to the exchange energy of the finite-range interaction. 

We recall here [221 EH] that without the tensor terms, relations f l6"i2"|) and f l6"3l allow 
us to determine the time-odd coupling constants Cf, C t As , and Cf as functions of the 
time-even coupling constants Cf, C^ p , and C t T . Since the time-even coupling constants 
are usually adjusted solely to the time-even observables, the resulting values of the 
time-odd coupling constants are simply "fictitious" or "illusory", as noted already in 
Ref. [3D]. In a more realistic case of relations ( 14"9"|) and fl5"D|) . these constraints are no 
longer valid, and the time-odd properties of the functional are independent of the time- 
even properties. This independence requires breaking the link between the Skyrme force 
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Table 1. The NV coupling constants (j49 |) -([50 | calculated for the Gogny interaction 
D1S [17] for the Fermi momenta of kp = and 1.35 fin -1 . First-order coupling 
constants, Cf and C t As , are in units of MeVfm 3 and the second-order coupling 
constants, C t Ap , C t r , C t As , and Cf, are in units of MeVfm 5 . 





t = 





t = 


1 




k F = 


k F = 1.35 


k F = 


kp = 1.35 




-665.1658 


-600.6156 


468.5360 


428.3580 


c; 


-25.09219 


-57.13246 


221.7219 


230.3318 


c? p 


-125.3365 


-84.66327 


56.65570 


31.22042 


c; 


236.9227 


74.22964 


-141.9368 


-40.19573 




10.72944 


-10.25276 


58.80425 


65.14281 


cj 


-80.70182 


3.226982 


-10.87263 


-36.22687 



Table 2. The standard Skyrme-force parameters (|5*T]) -([5l? ]) calculated for the Gogny 
interaction D1S [T7] for the Fermi momenta of kp — and 1 .35 fm -1 . These parameters 
correspond to the time-even sector of the Skyrme functional. Parameter to is in units 
of MeVfm 3 ; parameters t\ and £2 are in units of MeVfm 5 , and parameters xq, x\, and 
X2 are dimensionless. 





k F = 


k F = 1.35 


to 


-1773.775 


-1601.642 


h 


984.3584 


550.5103 


t 2 


810.3964 


-166.0710 


X 


0.5565848 


0.5697973 


Xi 


0.2488322 


0.09972511 


x 2 


-0.9915809 


-0.5517191 



and the density functional. 

4. Application to the Gogny interaction 

In this section, we apply the results of Sec. E] to the finite-range part of the Gogny 
interaction D1S [17]. This amounts to calculating moments (|52|) and (|53|) of the Gaussian 
functions with the two ranges of 0.7 and 1.2 fm, which constitute the central part of the 
Gogny interaction. Because this interaction does not contain any finite-range spin-orbit 
or tensor force in Eqs. ( 150]) and ( 151T) . the moments P, Q, R, and S are set to zero. On 
the other hand, the zero-range spin-orbit and density- dependent terms of the Gogny 
interaction are left unchanged. 

In Tabled] we show values of coupling constants (14"9]) - (I5"T]) calculated in the vacuum 
(kp = 0) and at the saturation density (kp = 1.35 fm -1 ). Similarly, Table |2] shows 
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Figure 2. Dependence of the NV coupling constants (|49|) - ([50l) . calculated for the 
Gogny interaction D1S |17) . on the Fermi momentum kp. Full and open symbols 
(lower and upper panels) show values of the time-even and time-odd coupling constants, 
respectively. Solid and dashed lines (left and right panels) show values of the isoscalar 
and isovector coupling constants, respectively. Units are specified in the caption to 
Table [0 



values of the Skyrme-force parameters (|57l) -(l59l ) corresponding to the time-even coupling 
constants. In Figs. [2] and |3] we plot the coupling constants and Skyrme-force parameters, 
respectively, as functions of the Fermi momentum. 

The most important observation resulting from values shown in Tables [1] and [2] 
and Figs. |2] and |3] pertains to a significant density (or kp) dependence of the coupling 
constants and Skyrme-force parameters. The strongest dependence is obtained for 
the isoscalar tensor coupling constant Cq . (Note that the central finite-range Gogny 
interaction induces significant values of the tensor coupling constants C f T , even if this 
force does not contain any explicit tensor term.) Also the kinetic coupling constants 
Cl exhibit a strong density dependence, going almost to zero at kp ~ 2fm _1 . We 
note that the obtained density dependencies do not, in general, follow any power 
laws. Significantly stronger density dependencies are obtained for the Skyrme-force 
parameters (Fig. [3]). The pole appearing in the parameter x 2 is a consequence of the 
fact that parameters t 2 and t 2 x 2 , derived in Eq. (|59|) . change signs at slightly different 
values of the Fermi momentum. 

Let us now consider the question of whether parameters calculated at any fixed 
value of the Fermi momentum can provide a reasonable alternative. To analyze this 
point, we first note that the zero-order coupling constants given in Eq. (l4"9l) depend on 
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Figure 3. Dependence of the standard Skyrme-force parameters (f57|) -(f59 |) . calculated 
for the Gogny interaction D1S [17], on the Fermi momentum kp. These parameters 
correspond to the time-even sector of the Skyrme functional. Units are specified in the 
caption of Table [2] 



kp (i) implicitly through the moments of Eq. ( |53|) and (ii) explicitly through the k 2 F 
term, that is, 

C?(k F ) = C t p t0 (k F ) + C7£ 2/3 (MP 2/3 , (65) 

c s t {k F ) = c s tfl {k F ) + q 2/3 (Mp 2/3 , (66) 

where we employed the standard association of the Fermi energy with density, namely, 
p = 2k F /37T 2 . Therefore, we can consider fixed values of the coupling constants 
calculated at a given Fermi momentum k F as 

C?(k° F ) and C s t (k° F ), or (67) 
Ct,o{k F ), C^2/ 3 (k F ), C^ Q (k F ), and C* )2 / 3 (/cp). (68) 

The second option gives the coupling constants that still depend on the density as p 2//3 , in 
analogy with the standard density-dependent term, which for the Gogny force depends 
on the density as p 1 ^ 3 , and which in the present study is alway kept untouched. 

In the upper and lower panels of Fig. HJ thin lines show the nuclear matter equations 
of state (energy per particle in function of density or Fermi momentum) obtained for 
the coupling constants fixed according to prescriptions ( IBTl) and (1681) . respectively. For 
comparison, thick lines show the NV results obtained for density-dependent coupling 
constants. Since in the nuclear matter, the factor multiplying the term z/ (r)^2(^) in 
Eq. ( Tl8l) vanishes exactly (by construction), the thick lines correspond to the exact 
Gogny force results. 

Equations of state calculated with prescription ( |68l) for fixed values of k° F from 1.35 
to 1.98 fm -1 with the step of 0.07 fm" 1 completely miss the saturation point. This means 
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Figure 4. Dependence of the infinite-matter energy per particle on the density 
(left panels) or Fermi momentum (right panels). Thick lines show the NV results 
for density-dependent coupling constants. Thin lines show the NV results for fixed 
coupling constants calculated at the given values of kP F . Upper and lower panels show 
results for fixed coupling constants (|67|) and for coupling constants (|68|) that depend 
on density as p 2 ^ 3 . The full square and circles mark the minima of curves. 



that Skyrme forces with constant parameters (l68j) derived by the NV expansion cannot 
be equivalent to the finite-range Gogny interaction. On the other hand, prescription 
(|67|) . for fixed value of k F = 1.35 fm -1 , reproduces the equation of state fairly well, with 
some deviations seen only at densities beyond the saturation point. For comparison, we 
also show results obtained with k F = 1 and 1.7fm _1 , which fit the equation of state at 
low and high densities, respectively. These results show that the explicit dependence 
on the Fermi momentum, which appears in the NV expansion, cannot be used to define 
the density dependence of the coupling constants. 

The full density dependence of all Skyrme-force parameters was recently 
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Figure 5. Similar to Fig. [2] except for the time-odd NV coupling constants only. 
The open symbols show exact results (pti?l) - (p)0|) as plotted in Fig. [21 whereas the full 
symbols show results inferred from the time-even sector by using Eqs. ([62]) and ((63]) ■ 
Solid and dashed lines (left and right panels) show values of the isoscalar and isovector 
coupling constants, respectively. 

Table 3. Similar to Table [T] except for the time-odd NV coupling constants inferred 
from the time-even sector by using Eqs. (|62l) and (j63[) for a Fermi momentum of 
1.35 fm" 1 . 





t = t=l 


c? s 
cj 


-27.94756 200.2052 
20.92673 23.21032 
-26.47083 -44.78631 



implemented in a spherical self-consistent code |31j . Here we use this implementation 
to test the results of the NV expansion against the full-fledged solutions known for 
the Gogny D1S force [32j [33] . In Table H] we show results obtained for three sets of 
Skyrme-force parameters: 

• SISa: fixed parameters given in Table |2]for k° F = 1.35 fm -1 , used together with the 
standard D1S parameters [17] of the zero-range spin-orbit (W = 130MeVfm 5 ) and 
density-dependent (£3 = 1390.60 MeVfm 4 ) terms. 

• SISb: density-dependent parameters shown in Fig. [3] implemented for p = 2kp/3Tr 2 
and used with the standard values of W and £3. 

• SISc: fixed parameters identical to those of SISa, but with the value of £3 = 
1385.35 MeVfm 4 slightly changed, so as to bring the ground-state energy of 208 Pb 
exactly to the true D1S value [33J. 

Table H] shows the ground-state energies E of seven doubly magic spherical nuclei, 
calculated by using the Skyrme-force parameters SISa, SISb, and SISc, and compared 
with the Gogny- force energies Eg- To facilitate the comparison, we also show relative 
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differences AE = (E — Eq)/\Eq\ in percent along with the corresponding RMS 
deviations (the last row in Table H]) . 

On can see that for the Skyrme-force parameters SISa and SISb, which are directly 
derived from the Gogny force by using the NV expansion, one obtains the nuclear binding 
energies smaller by 1-2% as compared to those given by the original Gogny force. This 
should be considered a very good result, although it cannot compete in precision of 
describing experimental data with the original Gogny or Skyrme forces, which have 
parameters directly fitted to experimental binding energies. A simple rescaling of the 
parameter i 3 brings the RMS deviation to below 0.4%, and makes the Skyrme force 
SISc competitive with most other standard Skyrme parameterizations. At the same 
time, the analogous RMS deviations obtained for the neutron and proton radii are 0.20 
and 0.30% (SISa), 1.01 and 0.94% (SISb), and 0.26 and 0.44% (SISc), respectively. 

We note here in passing that in this study we defined the SISb parameter set 
by considering the density- dependent coupling constants C t Ap (Fig. [3]) that multiply 
densities ptAp t . An attempt of using the same coupling constants along with densities 
— (Vp t ) 2 gives, in fact, the RMS deviations of binding energies (5.34%), neutron radii 
(2.44%), and proton radii (2.27%), which are significantly worse than those of SISb 
(Table HJ). This shows that the prescription to replace in EDFs integrated by parts (see 
Refs. [29| [31] ) the Fermi momentum by p = 2k F /3ft 2 may lead to significantly different 
results. Another recipy is to associate kp with density before taking products of the two 
density matrices which means both of the above terms will be active. However in order 
to have a better correspondence with the traditional Skyrme functionals we have made 
this association after taking the product. 

It is not our purpose here to propose that any of the Skyrme-force parameterizations 
introduced in the present work are better solutions to the problem of finding the 
best agreement with data. It is already known that within the standard second-order 
Skyrme-force parameterizations, a spectroscopic-quality [2] force cannot be found |34j. 
Nevertheless, it is gratifying to see that the NV expansion allows us to bridge the gap 
between the non-local and quasi-local EDFs, or between the finite-range and zero-range 
effective forces. A more quantitative discussion of the accuracy of the NV expansion 
will be possible by considering higher-order NV expansions [35] . 

Finally, in Fig. El we compare the time-odd coupling constants calculated by using 
Eqs. (j4"9~j) and (|50|) with those corresponding to the Skyrme-force parameters; that is, 
calculated by using Eqs. ( 1621) and ( 1631 . Similarly, Table [3] lists the numerical values of 
the Skyrme-force time-odd coupling constants. As one can see, differences between both 
sets of the time-odd coupling constants, shown in Fig. [5] with open and full symbols, are 
quite substantial. These results illustrate the fact that the NV expansion of the Gogny 
force leads to the Skyrme functional and not to the Skyrme force. 

We conclude this section by noting that functions fti{r) approximated by Gaussians 
(see Eqs. fl36|) - fl38|) and Fig. [I]) lead to the coupling constants and Skyrme-force 
parameters, which, when plotted in the scales of Figs. El EJ and El are indistinguishable 
from those presented in these figures. 
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Table 4. Binding energies E of seven doubly magic nuclei calculated by using the 
Skyrme-force parameters SISa, SISb, and SISc (see text) compared with the Gogny- 
force energies Eq. All energies are in MeV. 





D1S [33] 


SISa 


SISb 


SISc 


Eg 


E AE 


E AE 


E AE 


40 Ca 
48 Ca 
56 Ni 

100 Sn 
132 Sn 

208p b 


— 342.689 
— 414.330 
-481.111 
-637.845 
-828.024 
-1101.670 
-1637.291 


-335.312 2.15% 
-409.118 1.26% 
-473.497 1.58% 
-630.447 1.16% 
-814.568 1.63% 
-1086.272 1.40% 
-1612.634 1.51% 


-340.642 0.60% 
-410.698 0.88% 
-471.970 1.90% 
-629.066 1.38% 
-814.896 1.59% 
-1086.867 1.34% 
-1617.419 1.21% 


-339.369 0.97% 
-414.213 0.03% 
-479.843 0.26% 
-638.837 -0.16% 
-826.453 0.19% 
-1101.445 0.02% 
-1637.291 0.00% 


RMS 


n.a. 


n.a. 1.56% 


n.a. 1.33% 


n.a. 0.39% 



5. Conclusions 

In the present study, we derived a set of compact expressions giving the local energy- 
density corresponding to the Kohn-Sham potential energy for an arbitrary local finite- 
range central, spin-orbit, and tensor interactions. The method is based on the Negele- 
Vautherin density matrix expansion augmented by the odd-power gradient terms 
fulfilling the gauge-invariance condition. The coupling constants of the local energy 
density depend on a set of moments of the interaction conforming to the ideas of 
the effective theory. The expansion is based on the separation of scales between 
the range of the force and space characteristics of the one-body density matrix. It 
leads to a representation of dynamical properties of the system in terms of a set 
of numbers, whereby complicated short-range characteristics of effective interactions 
remain unresolved. 

We pointed out the fact that to correctly describe the exchange properties of 
the functional, proper treatment of the density matrix in the nonlocal direction is 
essential. This immediately leads to the local energy density that does not correspond 
to an averaged zero-range pseudopotential. Therefore, the Negele-Vautherin expansion 
performed up to NLO leads to the Skyrme functional and not to the Skyrme force. 
Within this formalism, the only way to define the Skyrme force is to match it to the 
time-even properties of the non-local functional, disregarding those pertaining to the 
time-odd channel. 

We applied the general NLO expressions to the case of the central finite-range part 
of the Gogny interaction. It turns out that the obtained coupling constants of the local 
Skyrme functional quite strongly depend on the Fermi momentum or on the density. 
Nevertheless, the equation of state obtained for fixed coupling constants, calculated 
at the saturation point of k,F = 1.35 fin -1 , fairly well reproduces the exact Gogny- 
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force result. On the other hand, partial density dependence, inferred from the explicit 
dependence of the coupling constants on hp, gives very unsatisfactory results. 

By solving the self-consistent equations with the Skyrme-force parameters derived 
from the Gogny force, one obtains an excellent agreement (up to 1-2%) of binding 
energies and radii with those corresponding to the true Gogny force. This shows that 
the ideas of the effective theory, whereby the finite-range nuclear forces are sufficiently 
short-range to be replaced by contact quasi-potentials, are applicable to low-energy 
nuclear observables. 

We also discussed properties of the one-body density matrix for the spin and isospin 
polarized infinite nuclear matter. In this case, one obtains the nonlocal densities with 
mixed spin-isospin channels. As a result, the local energy density is not invariant but 
covariant with respect to rotational and isospin symmetries; that is, it does not have the 
form of the standard Skyrme functional. It means that the standard Negele- Vautherin 
expansion can only be performed for unpolarized densities. 
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Contract No. N N 202 328234, by the Academy of Finland and University of Jyvaskyla 
within the FIDIPRO program, and by the U.S. Department of Energy under Contract 
Nos. DE-FC02-09ER41583 (UNEDF SciDAC Collaboration) and DE-FG02-96ER40963 
(University of Tennessee). 

Appendix A. Spin and isospin polarized infinite nuclear matter 

Let us consider the infinite nuclear matter with spin and isospin polarizations, which is 
described by the one-body Hamiltonian, 

h 2 

H = A-H<t-\ot, (A.l) 

2m K J 

where cr and H are the space vectors of spin Pauli matrices and spin-polarization 
Lagrange multipliers, respectively, and r and A are the isovectors of the analogous isospin 
Pauli matrices and isospin-polarization Lagrange multipliers, whereas the dot "•" (circle 
"o") denotes the scalar (isoscalar) product. Each eigenstate of Hamiltonian (lA.ljl is a 
Slater determinant that depends on the orientations of the Lagrange multipliers H and 
A in the space and isospace, respectively. However, since the kinetic energy is scalar 
and isoscalar, we can arbitrarily fix these orientations to H z and A3, which gives the 
nonlocal densities for spin-up and spin-down (cr = ±1) neutrons and protons (r = ±1) 
in the form, 

PaAr 1 ,r 2 )= f d 3 ke ik < r ^ =p ffT 7f 0T (r). (A.2) 

The system simply separates into four independent Fermi spheres for spin-up and spin- 
down neutrons and protons, with four constant densities p aT , whereas the dependence on 
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the relative position vector r is given by four scalar functions vr (JT (r) [compare Eq. I l30p ]. 

We note that the spin-isospin indices ar pertain here to the preselected quantization 
axis defined by the chosen directions of H and A, which define an "intrinsic" reference 
frame. In this reference frame, densities are marked with a bar symbol. 

The ground state of the system is obtained by filling the four Fermi spheres up to 
the common Fermi energy ep, 



2 U2 



h z k 



€p — 6p ( 



2m 



H o - A 3 r, (A.4) 
Finally, by varying ep, one obtains systems 



which defines the four Fermi momenta kp, 
with different total densities p = J^arPa-r- 

It is, of course, clear that for the asymmetric and polarized infinite nuclear matter, 
the density matrix of Eq. (1391) is diagonal in spin and isospin, 

pinatTi, r 2 a 2 r 2 ) = p aiT1 (r 1 , r 2 )5 aia2 S T1T2 (A. 5) 

and thus the nonlocal densities p^ifi, r 2) have non-zero components only for p, = or 
z and k = or 3, that is, 



( Poo{ri,r 2 ) \ 
Po3(r u r 2 ) 
pzo(ri,r 2 ) 
V Pz3(r h r 2 ) J 



( 1 
1 
1 

V i 



1 \ 

1 
1 

1 / 



( p ++ ( 

p+- 
p- 
V P- 



ri,r 2 ) \ 
ri,r 2 ) 
ri,r 2 ) 
ri,r 2 ) ) 



(A.6) 



pm{rx,r 2 ) 
Pzo(ri,r 2 ) 
V Pz 3 {ri,r 2 ) j 



\ **a(r) J 



(A.7) 



where we have abbreviated the indices of ar just to their signs. After expressing the 
right-hand side of this equation in terms nonlocal densities 0A.2p . one obtains: 

( Poo(ri,r 2 ) \ f 7f 00 (r) 7f 03 (r) Tt z0 {r) n z3 (r) \ ( p o \ 

*03(^) ttoo(^) Kzs{r) 7C z0 (r) po3 

K z o(r) n z3 (r) 7f 00 (r) 7f 03 (r) p z0 

V 7f*3(r) 7T 20 (r) 7r 03 (r) 7r 00 (r) / \ p z3 ) 
where functions vr At fc( r ) are defined similarly as in Eq. ( lA.6p . namely, 

/ ^oo(r) \ / 1 1 1 1 W vr ++ (r) \ 
1 1-1 1-1 vr + _(r) 
" 4 1 l_l_l vr_ + (r) 

V 1 -1 -1 1 A *—(r) ) 
Already here we see the main problem: for the spin and isospin polarized systems, 
the spin-isospin channels of nonlocal densities Pfj,k(vi,r 2 ) in Eq. (1 A. 7ft are linear 
combinations of the spin-isospin channels of local densities p^\ that is, the spin-isospin 
channels become mixed. 

To make the preceding result even more clear, we note that the spin-isospin 
directions of the Lagrange multipliers H and A can be arbitrarily varied and the spin- 
isospin directions of the nonlocal densities p^-fri, r 2 ), local densities p^, and functions 



(A.8) 
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^fc( r ) are always aligned with those of the Lagrange multipliers. Therefore, we can 
use the directions of the local densities instead of those pertaining to the Lagrange 
multipliers. By using the standard densities [19] in the (i) scalar-isoscalar channel 
(p = poo)) (h) vector-isoscalar channel (s M = p^o, for p = x,y,z), (hi) scalar-isovector 
channel (p& = pofc, for k = 1,2,3), and (iv) vector-isovector channel (s M & = p^t, for 
p = x,y,z and k — 1, 2, 3), we then define functions 7r(r) in the four channels as, 

vr(r) = 7f 00 (r), 
7r(r) = ri^Q^ir) , 

- g- (A-9) 

7?(r) = |vr z3 (r). 

Here, the "intrinsic" functions ^^(r) do not depend on the spin-isospin directions; that 
is, they are defined by the following Fermi energies, 

e F = e F>ar = _ \ H \ a _ |X| r . (A.10) 

Finally, definitions ( 1A.9|) allow us to present densities in the "laboratory" reference 
frame as [compare Eq. flA.7[) ] . 

p(fi, r 2 ) = pn{r) + p o 7f(r) + s ■ 7r(r) + s ■ 7r(r), ( A -H) 

^2) = p7r(r) + p7r(r) + s ■ 7r(r) + s • 7r(r), (A.12) 

s ( r i) r 2) = * 7r(r) + s o 7f(r) + pTt(r) + p o 7r(r), (A. 13) 

T2) = s 7r(r) + s 7r(r) + p7r(r) + p 7? (r). (A. 14) 

Note that the same scalar-isoscalar function 7r(r) multiplies all local densities in the 
first terms of Eqs. ( 1A.lll) -( lA.14l) . Therefore, the postulate of using different functions 
in different channels [11] is not compatible with the results obtained for the polarized 
nuclear matter. 

Again we see that the spin-isospin channels of nonlocal densities are mixed, namely, 
local densities in all channels contribute to every channel in the nonlocal density. As a 
consequence, the energy density is not invariant but only covariant with respect to the 
spin-isospin rotations (see the discussion in the Appendix A of Ref. [5]). Therefore, the 
NV expansion performed in the polarized nuclear matter does not lead to the standard 
local functional of Eq. ( )48l) . On the other hand, derivation in the unpolarized nuclear 
matter corresponds to all functions p aT { r ) equal to one another, which leads to vanishing 
functions vf(r), 7r(r), and n(r). Then, in Eqs. ( 1A.11I) -( |A.14I) . only the first terms survive 



and the spin-isospin channels are not mixed. Such a situation corresponds to postulating 
a channel-independent function 7r(r), which we employed in Sec. El 
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